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Abstract 

Background: The study describes the relationship of retinal vascular geometry (RVG) to severity of diabetic 
retinopathy (DR), and its predictive role for subsequent development of proliferative diabetic retinopathy (PDR). 

Methods: The research project comprises of two stages. Firstly, a comparative study of diabetic patients with 
different grades of DR. (No DR: Minimal non-proliferative DR: Severe non-proliferative DR: PDR) (10:10: 12: 19). 
Analysed RVG features including vascular widths and branching angles were compared between patient cohorts. 
A preliminary statistical model for determination of the retinopathy grade of patients, using these features, is 
presented. Secondly, in a longitudinal predictive study, RVG features were analysed for diabetic patients with 
progressive DR over 7 years. RVG at baseline was examined to determine risk for subsequent PDR development. 

Results: In the comparative study, increased DR severity was associated with gradual vascular dilatation (p = 0.000), 
and widening of the bifurcating angle (p = 0.000) with increase in smaller-child-vessel branching angle (p = 0.027). 
Type 2 diabetes and increased diabetes duration were associated with increased vascular width (p = <0.05 In the 
predictive study, at baseline, reduced small-child vascular width (OR = 0.73 (95% CI 0.58-0.92)), was predictive of 
future progression to PDR. 

Conclusions: The study findings suggest that RVG alterations can act as novel markers indicative of progression of 
DR severity and establishment of PDR. RVG may also have a potential predictive role in determining the risk of 
future retinopathy progression. 
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Background 

Diabetic retinopathy (DR) is a leading cause of visual 
loss in many developed countries [1,2]. Given that eyes 
with severe non proliferative diabetic retinopathy 
(NPDR) have a 52% risk of developing PDR within one 
year and 60% risk of developing high risk PDR within 
5 years [3] there is a need in clinical practice for effective 
risk stratification. Structural and functional changes in 
the retinal vasculature have been shown in several studies 
to be closely related to diabetes and DR [4,5] . Population- 
based cohort studies have described various vascular 
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calibre changes occurring with the development of 
diabetes [6], or in diabetics with no retinopathy [7], as 
well as with the development of early retinopathy in 
type 1 and type 2 diabetics [1,2,8,9], and indeed the 
progression of DR [7,10-13]. Some of the retinal vascular 
changes were also found to be significantly associated with 
progression risk to PDR while controlling for other 
baseline risk factors [14]. 

Changes in the retinal vascular calibre represent only a 
sole parameter of the retinal vascular network and do 
not convey information regarding the complexity of the 
retinal vascular branching pattern. According to Murray 
[15], the optimal vascular architecture achieves the most 
efficient blood flow transport with minimum energy 
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allowing for maximum vascular diffusion into the sur- 
rounding tissues. Alterations in the geometry of the retinal 
vascular network may thus reflect a state of vascular dys- 
function, and might potentially predict disease develop- 
ment [16-19]. Several features representative of the retinal 
microvasculature geometry have been evaluated in associ- 
ation with various systemic conditions [16-18,20-22]. 
Being non-dimensional, these features are less likely to be 
affected by variations in digital image resolution, ocular 
magnification or differences in refractive errors. In dia- 
betes, recent studies have explored the role of changes in 
retinal vascular geometry (RVG) as risk markers for inci- 
dent DR in young type 1 diabetics. Sasongko et al. [23] 
demonstrated that certain demographic factors such as 
age and sex as well as diabetes-related factors such as 
type, and the duration of diabetes were associated with 
subtle alterations in the RVG parameters including 
the branching angles. Moreover, other RVG features 
predicted incident retinopathy after a median of 3.8 years 
follow up [24] . 

In this study, we hypothesise that RVG changes are asso- 
ciated with increased severity of DR, and that, alterations 
in RVG can act as a novel marker in identifying networks 
at risk of future progression to proliferative retinopathy. 
We measure RVG at vascular bifurcations and analyse the 
relationship of these to progression of DR. 

Methods 

Study population 

The study was designed in two parts. Firstly a compara- 
tive observational study to evaluate the associations of 
RVG changes in age-matched cohorts of diabetic sub- 
jects with different grades of DR. Within this study, the 
relationships of certain demographic and clinical risk 
factors with RVG changes in the diabetic population 
were analyzed. Secondly, a pilot study of longitudinal 
data to assess the predictive role of RVG changes for 
subsequent development of proliferative retinopathy, 
using retrospective data for patients newly presenting 
with PDR. 

The studies were carried in accordance with the decla- 
ration of Helsinki (1989) of the World Medical Association. 
The protocols were approved by the local Sunderland 
Research Committee (SLREC 1129) and all recruited 
subjects gave written informed consent. 

Study participants were recruited from patients attend- 
ing the DR clinic at Sunderland Eye Infirmary and the 
local diabetic retinopathy screening programme centre. 

Type I or II diabetic patients aged 25 - 65 years, with 
different grades of DR and no previous pan retinal laser 
photocoagulation treatment were invited to participate in 
the study. The retinopathy grade was classified into; clinic- 
ally non-detectable, minimal non-proliferative, severe 
non-proliferative and proliferative diabetic retinopathy. 



Exclusion criteria included refractive error beyond +/- 
3 dioptres, concurrent ocular pathology including cata- 
ract, glaucoma, corneal opacities or any other retinal 
and optic nerve head pathology. Patients with uncon- 
trolled systemic hypertension with documented blood 
pressure measurements of more than 140/90 repeatedly 
recorded within the previous 12 months were also 
excluded. 

For the predictive study, patients who were referred 
with PDR from the diabetic retinopathy screening ser- 
vice to the hospital between January 2008 - January 
2009 were identified. Amongst these cases, patients 
aged 65 years or less who presented initially to the 
local diabetic retinopathy screening service with no 
clinically detectable retinopathy and with at least pre- 
ceding 6-years screening follow-up period were in- 
cluded (The progressors group). The control group for 
this study comprised of age matched diabetic patients 
attending the local screening service with no clinically 
detectable DR, also with a corresponding preceding 
6-year follow-up period and showing no signs of 
retinopathy progression through this period. (The non- 
progressors group). 

The other inclusion and exclusion criteria were similar 
to the comparative study. The baseline initial retinal 
screening images, the final screening images and the 
screening images captured at the penultimate visits were 
extracted for each eye and formally graded and analysed 
according to the study protocol. 

For all participants in both studies, one eye per sub- 
ject was randomly selected for analysis using online 
randomization programme. The demographic and medical 
data were extracted from the clinical notes. These included 
age, sex, type and duration of diabetes, history of hyperten- 
sion and high cholesterol as well as concurrent history of 
other diabetic microvascular complications. 

Photographic methods and retinopathy grading 

Mydriatic two-field 50° temporal and nasal retinal im- 
ages were captured for each eye by an experienced oph- 
thalmic medical photographer, according to the study 
protocol, using a Zeiss FF 450 Plus Fundus camera fitted 
with a JVC KY-F70B 3CCD digital camera with a reso- 
lution of 1360 x 1024 pixels. 

All images were assessed for combined field position 
and image quality according to the NHS diabetic eye 
screening programme image quality guidelines. Only im- 
ages with an overall "Good" quality scores were included 
for the study. 

DR was graded for these images by an experienced 
diabetic retinopathy grader, who was masked to the sub- 
jects' clinical data, according to the developed system for 
the EURODIAB IDDM complications study adaptation 
of the modified Airlie House classification [25]. 
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Retinal vascular geometry analysis 

The analysis of the bifurcation geometry was performed 
utilising a computer assisted semi-manual technique im- 
plemented on MATLAB software (MATLAB version 
7 0.1 Mathworks Inc. Massachusetts USA), as detailed 
elsewhere [26]. 

Briefly, for each retinal image, visible arteriolar and 
venular bifurcations were initially nominated by the ob- 
server across the whole image at varying distances from 
the optic disc, and above and below the horizontal mid- 
line. Care was also taken to include a widespread variety 
of bifurcations of different vascular size (Figure 1). Nom- 
inated bifurcations were then marked using the devel- 
oped tool; the vessel widths were measured by aligning 
rectangles over each vascular segment centreline and its 
width adapted to the observed segment edges. The junc- 
tion angles were estimated by inter-connecting four 
points; an intersection point and three end-points of the 
vessels centerlines (Figure 2). The tool presented the 
junction into a zoomed image allowing the measure- 
ments to be obtained at sub-pixel accuracy. For each 
image an average of 15 bifurcations were measured. 

For each analysed vascular bifurcation, the estimated 
structural and geometrical parameters were: the absolute 
width measurements of the parent, larger branch and 
smaller branch vessel segments (da di and <i 2 respectively); 
and the branching angles of the larger and smaller vessel 
segments, and the overall bifurcation angle (6x i 82 and 
6 = dj + 02 respectively) [16,27]. 

The precision of the custom-designed technique was 
compared to other semi-manual techniques used for 



vascular width detection, namely the "Kick points" [28], 
"Full width half maximum" [29], "Gregson's technique" 
[30], "Single Gaussian algorithm" [31], and the "Two 
dimensional Gaussian model" techniques [32]. The de- 
veloped rectangle technique' results were more precise 
than other techniques as reflected by lower standard de- 
viation and coefficient of variations values and tighter 
confidence intervals. Further details were presented in 
the author's thesis [33]. 

A single trained observer, masked to the subjects' 
demographic and clinical data, performed all the mea- 
surements in the study. The intra-observer repeatability 
was evaluated on a set of eight images (total of 117 
bifurcations) marked three times on separate occasions 
over a period of one month. The resultant Intra- 
observer repeatability was assessed using one-way ana- 
lysis of variance (ANOVA). The average coefficients of 
variation were 4.2% and 3.8% for the vascular width and 
angular measurements respectively. 

Statistical analysis 

Statistical analysis was performed using Minitab 16 stat- 
istical software. Chi Square test X 2 was used to detect 
differences between baseline characteristics. For the 
comparative study, one-way ANOVA was performed to 
compare between the subgroups of diabetic cohorts with 
different grades of diabetic retinopathy. A (p) value < 
0.05 was considered significant. In cases with significant 
(p) values, the Fisher-protected least significant differ- 
ence FLSD post-hoc test was used. The results were also 
sub analysed for arteriolar and venular subgroups. 
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Figure 2 Semi-manual technique for vascular analysis. The semi-manual technique: Vascular segment width estimated with aligned rectangles 
across the vessel longitudinal central axis. The bifurcation angles measured at the intersection point of vascular segments centrelines. 



Multiple linear regression analysis was then used to de- 
termine the effect of different demographic and clinical 
risk factors on the RVG parameters in the diabetic popula- 
tion with or without established retinopathy. Tested fea- 
tures included the continuous variables; age, duration of 
diabetes and categorical variables; sex, type of diabetes, 
history of hypertension, and history of hypercholesterol- 
emia. We constructed three models; model 1 included age 
and sex, model 2 additionally included the diabetes type 
and history of hypertension, model 3 additionally included 
diabetes duration and history of hypercholesterolemia. 

For the predictive study, one-way ANOVA was used to 
compare the RVG features between baseline, penultimate 
and final visits' images for both the progressors and non- 
progressors groups, and FLSD post-hoc test was per- 
formed for significant relations only where the (p) value 
was < 0.05. Binary logistic regression analysis was used to 
determine the odd ratios for retinopathy progression 
based on baseline RVG features measurements of both 
the progressors and non-progressors groups. 

For the diabetic retinopathy grade assessment, a leave- 
one-out cross validation creating logistic regression 
models for each patient. 



Results 

For the diabetic cross-sectional study, a total of 1518 
arteriolar and venular bifurcations of 51 eyes (51 subjects) 
were analysed. Details of the demographic and clinical 
baseline characteristics of the included subjects in the 
different diabetic subgroups are summarised in (Table 1). 

The groups were age-matched with no significant dif- 
ference in most of their demographic and clinical data. 
Unsurprisingly, the duration of diabetes was significantly 
longer in the PDR group as compared to the other 
subgroups, reflecting an established risk factor for 
progression of DR. 



Comparative study analysis 

For the comparative analysis between diabetic subjects 
with different DR grades, the distribution of measure- 
ments and estimates of the overall analysed retinal geo- 
metrical features and parameters in the different 
subgroups are presented in (Table 2). The results for the 
arteriolar and venular bifurcations subgroups are pre- 
sented in (Table 3). The overall significance of ClXthe 
ANOVA test is shown, along with the significant differ- 
ences in-between the subgroups as determined by the 
FLSD test. 

The results show that increased DR severity was asso- 
ciated with greater vascular segment width. These rela- 
tionships are maintained for the arteriolar and venular 
bifurcations subgroups (p = 0.000). 

There is also a widening of the total bifurcation angle 
(8) with increased severity of DR and development of 
proliferative retinopathy (p = 0.000). This widening 
reflected specifically an associated steady increase in the 
branching angle of the smaller-child vessel segment (d 2 ) 
with DR progression (p = 0.027). Such relationships for 
(6) and (8 2 ) angles were maintained for the arteriolar 
bifurcations (p = 0.000, p = 0.014 respectively); however 
in the venular bifurcations, the same trend of angular 
widening did not reach a statistically significant level 
(p = 0.062, p = 0.504 respectively). 

Relationship of RVG with demographic and clinic risk 
factors 

The relationship between the different demographic and 
clinical risk factors for disease progression and the vari- 
ous RVG features was demonstrated with the multiple 
regression analysis for the overall measurements and the 
sub-analysed arteriolar and venular subgroups. 

Increased age was associated with significant de- 
creased arteriolar and venular vascular width (p = <0.01) 
together with reduction in JER (p = <0.01). 
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Table 1 Demographic and clinical baseline characteristics of diabetic subjects 


Demographic and clinical baseline characteristics of diabetic subjects 




No retinopathy group 


Minimal NPDR group 


Severe NPDR group PDR group 


P value 


Subjects 


10 Subjects 


10 Subjects 


1 2 Subjects 


19 Subjects 




(Eyes) 


(10 eyes) 


(1 0 eyes) 


(12 eyes) 


(19 eyes) 




Age 
Mean 


55.5 


56.6 


53 


50 














0.92 








(9 01 


(12 5) 




(Range) 




p/-coj 


f41 - 651 


(26-65) 






60% 


50% 


42% 


58% 


021 


Type Diabetes 


70% 


80% 


83% 


73% 


0.70 


Duration of DM in months: Mean 


118 


148 


128 


211 




(SD) 


(65) 


(74) 


(93) 


(124) 


<0.05 


(Range) 


(24-144) 


(37-240) 


(1 2-300) 


(2-456) 




Hypertension 


50% 


70% 


67% 


48% 


0.74 


High Cholesterol 


40% 


70% 


34% 


57% 


0.16 


Associated microvascular complications 










0.44 


NONE 


80% 


70% 


59% 


63% 




IHD, CVA, TIA 


10% 


20% 


0% 


5% 




Neuropathy, Foot Ulcer 


10% 


10% 


25% 


10% 




Nephropathy 


0% 


0% 


17% 


21% 




Smoking 


20% 


30% 


17% 


21% 


0.8 


Duration of DM has a statistically significant correlation with Type of DR, as expected; the other factors do not. 






Female sex was associated with decreased venular 


vas- was associated with increased both the arteriolar bifur- 


cular width (p = <0.05) (apart from the diameter of the cation angle 


(d) (p = <0.05) 


and the smaller-child 


smaller-child vessel segment) as well as reduction in 


ar- branching ang 


;le (B 2 ) (p = <0.05). 






teriolar and venular JER (p = <0.05). 










Type 2 diabetes was associated with increased arteri- Model for diabetic retinopathy grade assessment 




olar vascular width (p = <0.05) while the increased dur- We conducted a preliminary study to determine if a stat- 


ation of diabetes was associated with increased arteriolar istical model 


can determine the grade of retinopathy 


and venular vascular width (p = 


= <0.01). 


given the RVG feature distribution for a particular 


image 


History of hypertension was 


associated with increased (and therefore patient). We used the means of the RVG 


venular vascular width (p = <0.05) (excluding the small- features to characterize each imai 


ge. Given the small size 


child vessel segment), while history of high cholesterol of the data set (51 patients), we 


used the leave-one-out 


Table 2 The distribution of geometrical measurements in the diabetic subgroups 






Overall data 


Retinal parameter 


No retinopathy 


Mild NPDR group 


Severe NPDR 


PDR group 




Mean ± (STD) 


group 


group 


P 

value 




(242 bifurcations) 


(310 bifurcations) (372 bifurcations) (594 bifurcations) 




Parent Vessel Diameter d 0 (pixels) 


6.97 ± (1.61)* 


8.39 ± (2.37) t 


8.56 ± (2.39)+ 


8.98 ± (2.63)+ 


0.000 


Large Child Diameter d, (pixels) 


6.17 ± (1.56)* 


7.44 ± (2.30)t 


7.56 ± (2.34)+ 


8.01 ± (2.52)+ 


0.000 


Smaller Child Diameter d 2 (pixels) 


4.66 ± (1.04)* 


5.47 ± (1.53)t 


5.51 ± (1.44)+ 


5.76 ± (1.70)+ 


0.000 


Bifurcating Angle 6 (Degrees) 


77.04 ± (15.6)* 


79.85 ± (17.4)*t 


80.02 ± (17.8)+ 


83.23 ± (18.9)+ 


0.000 


Branching Angle 0, (Degrees) 


25.14 ± (14.1) 


25.12 ± (15.8) 


25.82 ± (15.2) 


25.05 ± (16.7) 


0.941 


Branching Angle 9 2 (Degrees) 


51.97 ± (19.3)* 


55.01 ± (21.2)* 


54.84 ± (21.7)*+ 


58.57 ± (23.7)+ 


0.027 


Junction Exponent x 


3.33 ± (3.24)*+ 


3.34 ± (3.17)*+ 


3.30 ± (3.13)+ 


3.55 ± (3.32)* 


0.050 



The p-value of ANOVA test is shown with the significant differences between the subgroups as determined by the FLSD test. (Means that do not share a symbol 
are significantly different). 
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Table 3 The distribution of the arteriolar and venular geometrical measurements in the diabetic subgroups 



Arteriolar data 



Retinal parameter 


No retinopathy group 


Mild NPDR group 


Severe NPDR group 


PDR group 


p value 


Mean ± (STD) 


(117 bifurcations) 


(136 bifurcations) 


(139 bifurcations) 


(231 bifurcations) 




Parent Vessel Diameter d 0 (pixels) 


6.47 ± (1.08)* 


7.54 ± (1.68)t 


7.58 ± (1.65)t 


8.06 ± (1.84)* 


0.000 


Large Child Diameter d, (pixels) 


5.78 ± (1.10)* 


6.65 ± (1.62)t 


6.73 ± (1.66)t 


7.23 ± (1.80)* 


0.000 


Smaller Child Diameter d 2 (pixels) 


4.60 ± (0.99)* 


5.42 ± (1.37)t 


5.43 ± (1.33)t 


5.62 ± (1.43)t 


0.000 


Bifurcating Angle 6 (Degrees) 


76.13 ± (15.8)* 


79.14 ± (17.9)* 


78.78 ± (18.4)* 


84.33 ± (18.6)t 


0.000 


Branching Angle 0, (Degrees) 


25.77 ± (15.2) 


27.32 ± (16.6) 


28.33 ± (15.4) 


27.44 ± (16.4) 


0.560 


Branching Angle 9 2 (Degrees) 


50.36 ± (20.3)* 


51.99 ± (21.3)* 


50.45 ± (20.3)* 


56.88 ± (22.9)t 


0.014 


Junction Exponent x 


3.70 ± (1.27) 


3.80 ± (2.00) 


3.92 ± (1.85) 


4.15 ± (2.25) 


0.161 


Venular data 


Retinal parameter 


No retinopathy group 


Mild NPDR group 


Severe NPDR group 


PDR group 


p value 


Mean ± (STD) 


(125 bifurcations) 


(1 74 bifurcations) 


(233 bifurcations) 


(363 bifurcations) 




Parent Vessel Diameter d 0 (pixels) 


7.42 ± (1.87)* 


9.04 ± (2.62)t 


9.13 ± (2.56)+* 


9.55 ± (2.88)* 


0.000 


Large Child Diameter d, (pixels) 


6.51 ± (1.83)* 


8.04 ± (2.55)t* 


8.05 ± (2.54)t 


8.50 ± (2.79)* 


0.000 


Smaller Child Diameter d 2 (pixels) 


4.72 ± (1.08)* 


5.51 ± (1.64)t 


5.55 ± (1.47)t 


5.82 ± (1.85)* 


0.000 


Bifurcating Angle 6 (Degrees) 


77.68 ± (15.1) 


80.40 ± (16.9) 


80.75 ± (17.4) 


82.53 ± (19.1) 


0.062 


Branching Angle 0, (Degrees) 


24.81 ± (12.9) 


23.43 ± (14.9) 


24.34 ± (15.0) 


23.53 ± (16.8) 


0.921 


Branching Angle 6 2 (Degrees) 


53.00 ± (1 7.9) 


57.34 ± (20.9) 


57.43 ± (22.1) 


59.65 ± (23.7) 


0.504 


Junction Exponent x 


3.02 ± (0.96) 


2.99 ± (0.88) 


2.94 ± (0.99) 


3.18 ± (1.46) 


0.07 



The p-value of ANOVA test is shown with the significant differences between the subgroups as determined by the FLSD test. (Means that do not share a symbol 
are significantly different). 

There are significant differences in the vessel widths in both arteriolar and venular data, and in the second child branching angle, 8 2 , in arteriolar data. 



cross validation methodology, constructing logistic re- 
gression models to predict the grade of each patient, and 
training each model using the other patients. We evalu- 
ated multinomial models to predict all four grades, but 
found it more effective to consider two separate types of 
logistic regression model: model one distinguishes pa- 
tients without retinopathy from those with retinopathy, 
and model two distinguishes non-proliferative from pro- 
liferative retinopathy. 

The RVG parameters vary systematically with distance 
from the ONH. However, the distribution of bifurcation 
distances may vary according to the clarity of the images 
(with a larger proportion further from the ONH for 
clearer images), and the clarity in turn may be correlated 
with grade. To ensure that the experiments are con- 
trolled for distance from ONH we applied the following 
methodology. First, we identified the patient with the 
lowest mean distance ratio from ONH (i.e. distance over 
ONH diameter). Second, for each of the other patients 
we progressively discarded the most distant (by ratio) bi- 
furcations until the mean distance ratio was just lower 
than the "minimum" patient. This procedure ensured 
that the mean ONH distance ratios for patients were all 
approximately equal. 

Model one achieved 97.6% sensitivity and 90% specifi- 
city to determine the presence of retinopathy versus no 



retinopathy, using two predictors: mean parent diameter 
and mean deflection angle of the smaller child. Adding 
other features did not improve the model accuracy. Util- 
izing the mean parent diameter only reduced specificity 
to 80%, sensitivity remained at 97.6%, indicating that the 
mean deflection angle made a small contribution to 
model performance. Utilizing mean deflection angle 
alone did not yield a usable model. For this model we 
utilized all 51 patients, with the mild, severe and prolif- 
erative classes merged into a single "retinopathy" class. 

Model two achieved 63.2% sensitivity and 72.3% speci- 
ficity to determine proliferative versus non-proliferative 
retinopathy. This model used mean parent diameter as 
the sole predictor; adding the mean deflection angle did 
not improve performance. For this model we utilized the 
41 retinopathy cases, distinguishing the 22 non- 
proliferatives from 19 proliferatives. If the 10 "no retin- 
opathy" cases are also included in the study the model 
correctly classifies them as non-proliferative, yielding an 
overall specificity of 81.3%. 

Predictive study analysis 

For the predictive study, five patients aged 50 - 65 years 
(mean = 57 years) were identified with PDR who fulfilled 
the inclusion criteria. The non-progressors group 
included five patients aged 47 - 65 years (mean = 59 years). 
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A total of 121 bifurcations for the progressors and 97 
bifurcations for the non-progressors were analysed per 
screening visit. 

The distributions of the geometrical measurements 
and estimates in the progressors group for the three 
screening visits are shown in (Table 4) together with the 
results of the one-way ANOVA test. Absolute vascular 
width comparisons were excluded from the analysis due 
to significant difference in the used camera systems over 
the follow-up period with difference in image resolution 
that could hamper comparisons made in terms of abso- 
lute pixel size. The rest of the geometrical measure- 
ments, being non-dimensional, thus unaffected by image 
resolution, are presented. 

There was no statistically significant difference noted 
in any of the analysed RVG features in the progressive 
group between the three screening visits, however; there 
was a noted widening of the vascular bifurcating angle 
(6) together with a corresponding widening of the 
smaller-child vessel branching angle (6 2 ) in the final visit 
with development of PDR - a trend similar to that ob- 
served in the cross-sectional comparative study. Such 
trend was maintained for the arteriolar and venular sub- 
group analysis (Table 5). 

For the non-progressors group, no statistically signifi- 
cant difference was noted for any of the analysed RVG 
parameters between the three screening visits, and no 
obvious trend in angular width changes could be ob- 
served similar to the progressors group. 

The predictive value of RVG features at baseline for 
determining future progression to PDR was evaluated 
using the binary logistic regression analysis. The results 
are shown in (Table 6). Reduced arteriolar and venular 
smaller-child vessel width (d.2) at baseline significantly 
predicted future progression to PDR (OR = 0.73 (95% CI 
0.58-0.92)). The results reached statistical significance in 
most of the subgroup analysis. 

Discussion 

In this study, we evaluated the association of RVG alter- 
ations with increased severity of DR in age-matched co- 
horts of adult type 1 and 2 diabetics with established DR 
ranging from minimal NPDR to PDR. 



The main findings of our comparative study revealed a 
significant increase in all of the arteriolar and venular 
vascular segments width measurements with increased 
severity of DR. This was associated with a larger mean 
deflection angle of the smaller child vessel at bifurca- 
tions. The association of increased vascular calibre with 
retinopathy severity is consistent with results from previ- 
ous studies which demonstrated vasodilation in adults 
with type 2 diabetes with increasing diabetes duration. 
This was proposed to result from vessel wall smooth 
muscle cell and pericyte loss with vasomotor dysregula- 
tion and loss of structural support leading to the sus- 
tained dilatation [4,34,35]. The cause of the apparent 
change in mean bifurcation angle is not known, but may 
relate to either a departure from optimality or to adapta- 
tion to sustain optimality. These results are consistent 
with Sasongko et al. showing larger arteriolar branching 
angles in association with longer diabetes duration [23], 
but further suggest that this is primarily due to increased 
deflection of the smaller child vessel. 

The results from our patient level statistical model 
study must be treated with caution, given the small 
sample size. They are indicative that the distribution of 
vessel diameter measurements, taken at bifurcations, 
may be useful in identifying the presence of retinopathy 
and, to a lesser extent, proliferative retinopathy. The 
mean deflection angle of the child vessels may also be of 
some predictive value, although this is less clear. 

In the predictive study, we attempted to establish the 
clinical relevance of the comparative study findings on 
an individual level. The results showed the same trend, 
previously demonstrated in the comparative study, of bi- 
furcation angular widening with progression of retinop- 
athy and reaching the maximum with established PDR 
after an average of 6 years follow up. These changes were 
rather clearer in the arteriolar vascular system than the 
venular system, yet did not reach statistical significance, 
which might be related to the small number of recruited 
subjects. 

The predictive study results have also shown that in 
adult diabetics with initially no clinically detectable ret- 
inopathy, alterations in RVG at baseline in the form of 
narrower arteriolar and venular small-child vessel 



Table 4 Overall Geometrical features of progressors group 



Progressors overall data 


Retinal parameter 


Baseline visit 


Penultimate visit 


Final visit 


P 


Mean ± (STD) 


(No retinopathy) 


(No or minimal NPDR) 


(PDR) 


value 


Branching Angle 0, (Degrees) 


2647 ± (14.9) 


28.49 ± (15.3) 


26.26 ± (13.9) 


0.501 


Branching Angle 9 2 (Degrees) 


51.98 ± (20.5) 


49.98 ± (21.0) 


54.07 ± (17.6) 


0.334 


Junction Exponent x 


3.102 ± (1.02) 


3.01 ± (0.83) 


3.104 ± (1.08) 


0.721 



Means and Standard deviations for the geometrical features in the progressors group baseline, penultimate and final visits for the overall data. P value for ANOVA 
test is shown. 
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Table 5 Arteriolar and venular geometrical features for progressors group 



Progressors arteriolar data 


Retinal parameter 
Mean ± (STD) 


Baseline visit 


Penultimate visit 


Final visit 


p value 


Branching Angle 0, (Degrees) 


24.94 ± (14.0) 


25.77 ± (11.6) 


24.71 ± (11.9) 


0.924 


Branching Angle 6 2 (Degrees) 


54.58 ± (22.8) 


54.47 ± (22.1) 


57.87 ± (18.7) 


0.675 


Progressors venular data 


Retinal parameter 
Mean ± (STD) 


Baseline visit 


Penultimate visit 


Final visit 


p value 


Branching Angle 0, (Degrees) 


27.72 ± (15.7) 


27.72 ± (15.4) 


27.47 ± (15.2) 


0.995 


Branching Angle 9 2 (Degrees) 


49.87 ± (1 8.4) 


47.43 ± (17.7) 


51.10 ± (16.2) 


0.561 



Means and Standard deviations for the geometrical features in the progressors group baseline, penultimate and final visits for the arteriolar and venular data. 
P value for ANOVA test is shown. 



segment diameter at a bifurcation and a reduced asym- 
metry and area ratio values may be correlated with fu- 
ture progression to PDR. These results echo the earlier 
work by Fanucci et al. who suggested that architectural 
changes in the arterial bifurcations of human vascular 
networks in the form of narrow angles and normal area 
ratios or wide angles and reduced area ratios may pre- 
dispose to future vascular diseases [36]. 

Our results also demonstrate the effect of different 
demographic and clinical risk factors on the RVG fea- 
tures in adult diabetic subjects with various stages of 
existing retinopathy. Age was associated with reduced 
arteriolar and venular vascular widths, corresponding to 
earlier work by Klein et al. on type 1 and type 2 diabetics 
[7,37]. Longer duration of diabetes was associated with 
arteriolar and venular dilatation. The increased arteriolar 
bifurcating angles noted here with high cholesterol levels 
supports the evolving evidence in previous studies of the 
role of lipids on the retinal microvasculature [38]. 

In this study, type 1 and type 2 diabetic patients were 
included for analysis. Type 2 diabetes was found to be 
associated only with increased arteriolar vascular width, 
with no other bearing on the rest of the analysed RVG 
features. There are no other previous reports that sug- 
gested definite comparative difference in retinal vascular 



analysis between both types of diabetes. In fact, in a re- 
cent study evaluating neural and vascular alterations in 
adolescent type 1 and type 2 diabetic retinae, there was 
no statistically significant difference in the arteriolar and 
venular widths, retinal function as measured by multi- 
focal electroretinogram or retinal thickness as detected 
by the optical coherence tomography noted between 
both diabetic groups [39]. 

It is therefore worth noting here that the observed al- 
terations in RVG associated with the various risk factors 
would have not affected our previous comparative and 
predictive results, as the study diabetic subgroups were 
evenly matched in most of their demographic and 
clinical data (Table 1). 

To the best of our knowledge, this is the first study to 
report on the RVG alterations in adult diabetics with 
various degrees of DR from mild retinopathy to estab- 
lished proliferative retinopathy. The study provides 
cross-sectional comparative analysis as well as a pilot 
predictive analysis for these variations with increased 
retinopathic severity. Most previous studies have used 
other semi-manual techniques of retinal vascular mea- 
surement that have focused on limited zones surrounding 
the optic disc [40]. In contradistinction, this study ana- 
lysed vascular bifurcations of different size and orientation 



Table 6 Results of binary logistic regression analysis 



Binary logistic regression analysis 







Overall 






Arteriolar 






Venular 




Retinal parameter 


Coeff 


OR (95% CI) 


p value 


Coeff 


OR (95% CI) 


p value 


Coeff 


OR (95% CI) 


p value 


do 


-0.025 


0.97 (0.84 - 1.14) 


0.742 


-0.145 


0.86 (0.64 - 1.18) 


0.352 


-0.021 


0.98 (0.81 - 1.18) 


0.831 


d, 


-0.007 


0.99 (0.84 - 1.17) 


0.93 


-0.178 


0.84 (0.61 - 1.16) 


0.279 


0.031 


1 .03 (0.84 - 1 .27) 


0.768 


d 2 


-0.318 


0.73 (0.58 - 0.92) 


0.007 


-0.358 


0.70 (0.49 - 1 .00) 


0.05 


-0.301 


0.74 (0.55 - 1 .00) 


0.05 


e, 


-0.001 


1 .00 (0.98 - 1 .02) 


0.844 


-0.007 


0.99 (0.97 - 1 .02) 


0.574 


0.002 


1 .00 (0.98 - 1 .03) 


0.858 


e 2 


0.011 


1.01 (1.00- 1.02) 


0.154 


0.014 


1.01 (1.00 - 1.03) 


0.125 


0.005 


1.01 (0.98 - 1.03) 


0.632 


X 


-0.118 


0.83 (0.68- 1.01) 


0.069 


-0.181 


0.83 (0.66 - 1 .06) 


0.131 


-0.099 


0.90 (0.54 - 1 .53) 


0.709 



Results of binary logistic regression analysis for the overall data, arteriolar and venular subgroups. The coefficient, odds ratios, 95% confidence intervals, and p 
values are shown. (Significant results in Bold). 



Habib et al. BMC Ophthalmology 2014, 14:89 
http://www.biomedcentral.com/1471-2415/14/89 



Page 9 of 1 1 



across the entire retina that could arguably be more repre- 
sentative of the retinal vascular network state in normal 
and disease states, especially with peripheral retinal is- 
chaemia associated with progressive DR. The study's vas- 
cular measurements were performed by a robust and 
reliable technique and evaluations were conducted by a 
single observer, masked to the retinopathy grade, to 
eliminate the inter-observer variability of measurements 
and reduce sources for bias. 

The limitations of the study reside in the fact that, 
despite our efforts for accurate data collection, it relies 
on a limited amount of retrospective clinical data. The 
computer-assisted semi-manual technique was time con- 
suming and labour intensive and thus limited the number 
of subjects. Future larger-scale studies, utilising reliable 
and more robust fully automated techniques of vascular 
analysis, are needed which would allow detailed statistical 
analysis to be conducted on large samples at the 
bifurcation-level and the patient-level. The results of the 
pilot predictive study are also limited by the small number 
of included patients that fitted the inclusion criteria; how- 
ever this reflects the rarity of such patients as previously 
demonstrated in the EURODIAB complication study 
where only 4% of subjects without diabetic retinopathy at 
baseline were associated with progression to PDR after 
follow-up period of 7.3 years [41]. We also note that the 
expected RVG parameters distribution varies with dis- 
tance from the ONH, and this distribution may vary 
between retinopathy grades for reasons other than patho- 
logical changes to individual bifurcations, including image 
clarity. Similarly, measurements such as vessel widths may 
be affected by media opacity. Although we controlled for 
the distance from the ONH in our patient level models, 
the approach does not admit for absolute consistency 
between patients having different retinal vascular geometry, 
making it difficult to ensure that the results do not partially 
reflect some measurement and/or statistical artifacts. 

Finally, in our results, it was noted that JER values suf- 
fered from over-sensitivity and wide scatter. This over- 
sensitivity of its calculations to minor variation in vessel 
width estimations might explain our wide range of mea- 
surements and high variability and make it a less than 
ideal value for investigational studies. 

Conclusion 

In conclusion, the study presents novel retinal vascular 
geometrical markers indicative of progression of DR se- 
verity and the establishment of proliferative retinopathy. 
The results suggest that these markers might allow fu- 
ture quantitative assessment of diabetic microvascular 
damage and have a potential predictive role in determining 
the risk of future retinopathy progression. Such findings 
could help reclassifying patient risk, stratify and alter their 
follow-up and management plans, and promote early 



identification of patients at risk of visual loss. The findings 
warrant further prospective studies to confirm and sup- 
port these results, aided by the ongoing developments in 
fully automated vascular analysis software packages. 

Abbreviations 

ANOVA: Analysis of variance; DR: Diabetic retinopathy; ETDRS: Early treatment 
diabetic retinopathy study; FLSD: Fischer-protected least significant 
difference; JER: Junction exponent ratio; LDR: Length/Distance ratio; 
NPDR: Non-proliferative diabetic retinopathy; ONH: Optic nerve head; 
PDR: Proliferative diabetic retinopathy; RVG: Retinal vascular geometry; 
WESDR: Wisconsin epidemiological study of diabetic retinopathy. 

Competing interests 

The authors declare that they have no competing or financial interests. The 
authors received no external funding support for collecting and analysing 
the data for this study. 

Authors' contributions 

MH has been involved in designing the study and acquisition of the images 
datasets, semi-manual marking and analysis of the retinal vascular geomet- 
rical features, statistical analysis of data, and drafting the manuscript of this 
study. BA made extensive contributions in adjusting the custom-designed 
vascular analysis tool, extracting and interpreting the different vascular 
features and analysing the data. AH made substantial contribution in 
developing the computerised software tool, designing the predictive model, 
analysing the results and drafting and revising the manuscript. DS conceived 
the study design, provided scientific supervision for the analysis and 
interpretation of the results, together with extensive contribution to the 
manuscript draft and critical analysis. All authors read and approved the 
final manuscript. 



Authors' information 

Mr Maged Habib is a consultant Ophthalmologist and Vitreoretinal Surgeon 
at Sunderland Eye Infirmary in the North East of England. He has special 
clinical interest in medical and surgical retinal diseases, specially the 
management of diabetic eye disease and retinal vascular disorders. He is the 
lead for virtual diabetic retinopathy surveillance service at South of Tyne 
Diabetic Eye screening programme. His main research projects explored the 
role of retinal vascular geometry analysis in the detection and monitoring of 
diabetic retinopathy as well as other systemic and ocular diseases. 
Dr. Bashir Al-Diri is a senior lecturer in School of Computer Science at the 
University of Lincoln. His main research is concerned with the development 
of computer vision algorithms to analyse retinal blood vessels for 
characterization of diabetic retinopathy. Dr. Al-Diri developed a robust fully 
automated system for retinal vascular segmentation and measurement, 
which provides a unique combination of good segmentation and superior 
measurement performance. It is thus uniquely well-suited to act as a base 
for research into the diagnosis of vascular diseases that cause measurable 
changes to the geometry of retinal vessels. Dr. Al-Diri research includes 
computer vision, medical image analysis, automated surveillance, artificial 
intelligence, speech recognition, language corpus and lexical analysis, and 
applications of active contour models. 

Prof. Andrew Hunter is Pro Vice Chancellor, Head of the College of Science, 
at the University of Lincoln. He has published over 100 academic papers, 
including many in international journals. He has also developed several 
freeware and commercial artificial intelligence software packages. His major 
research interest is in medical applications of computer vision. He is the 
Coordinator for the EU Marie-Curie ITN Project "REVAMMAD," which aims to 
integrate techniques for retinal vascular modelling, measurement and 
diagnosis. 

Mr David Steel is a Consultant Ophthalmologist at Sunderland Eye Infirmary 
in the North East of England. He is a vitreoretinal surgeon with clinical 
interests in all aspects of retinal disease and complex cataract surgery. He is 
an honorary Senior lecturer at the Institute of Genetic Medicine at University 
of Newcastle upon Tyne. Mr Steel is also the Consultant lead for the south of 
the Tyne Diabetic Retinopathy Eye-screening Programme. 



Habib et al. BMC Ophthalmology 2014, 14:89 
http://www.biomedcentral.com/1471-2415/14/89 



Page 1 0 of 1 1 



Acknowledgments 

We would like to thank the medical team at the diabetic retinopathy clinics 
at Sunderland Eye Infirmary and the South of the Tyne diabetic retinopathy 
screening programme for providing the retinal images and assisting in 
collecting the medical data for the recruited subjects. 

Author details 

'Sunderland Eye Infirmary - Queen Alexandra Road, Sunderland SR2 9HP, UK. 
department of Computing and Informatics, Lincoln University, Brayford 
Pool, Lincoln LN6 7TS, UK. institute of Genetic Medicine, University of 
Newcastle Upon Tyne, Central parkway, Newcastle NE1 3BZ, UK. 

Received: 3 March 2014 Accepted: 30 June 2014 
Published: 7 July 2014 

References 

1 . Congdon NG, Friedman DS, Lietman T: Important causes of visual 
impairment in the world today. JAMA 2003, 290(15)2057-2060. 

2. Fong DS, Aiello L, Gardner TW, King GL, Blankenship G, Cavallerano JD: 
Diabetic retinopathy. Diabetes Care 2004, 27(1 0):2540-2553. 

3. Early photocoagulation for diabetic retinopathy. ETDRS report number 9. 
Early Treatment Diabetic Retinopathy Study Research Group. 
Ophthalmology 1991,98(5 Suppl):766-785. 

4. Curtis TM, Gardiner TA, Stitt AW: Microvascular lesions of diabetic 
retinopathy: clues towards understanding pathogenesis? Eye (Lond) 2009, 
23(7):1496-1508. 

5. Mandecka A, Dawczynski J, Vilser W, Blum M, Muller N, Kloos C, Wolf G, 
Muller UA: Abnormal retinal autoregulation is detected by provoked 
stimulation with flicker light in well-controlled patients with type 1 
diabetes without retinopathy. Diabetes Res din Pract 2009, 86(1)51-55. 

6. Wong TY, Klein R, Sharrett AR, Duncan BB, Couper DJ, Tielsch JM, Klein BE, 
Hubbard LD: Retinal arteriolar narrowing and risk of coronary heart 
disease in men and women. The Atherosclerosis Risk in Communities 
Study. JAMA 2002, 287(9): 1 153-1 159. 

7. Klein R, Klein BE, Moss SE, Wong TY, Sharrett AR: Retinal vascular caliber in 
persons with type 2 diabetes: the Wisconsin Epidemiological Study of 
Diabetic Retinopathy: XX. Ophthalmology 2006, 113(9):1 488-1 498. 

8. Alibrahim E, Donaghue KC, Roger S, Hing S, Jenkins AJ, Chan A, Wong TY: 
Retinal vascular caliber and risk of retinopathy in young patients with 
type 1 diabetes. Ophthalmology 2006, 113(9):1499-1503. 

9. Gardiner TA, Archer DB, Curtis TM, Stitt AW: Arteriolar involvement in the 
microvascular lesions of diabetic retinopathy: implications for 
pathogenesis. Microcirculation 2007, 14(1)25-38. 

10. Islam FM, Nguyen TT, Wang JJ, Tai ES, Shankar A, Saw SM, Aung T, Lim SC, 
Mitchell P, Wong TY: Quantitative retinal vascular calibre changes in 
diabetes and retinopathy: the Singapore Malay eye study. Eye (Lond) 
2009, 23(8):1719-1724. 

11. Klein R, Klein BE, Moss SE, Wong TY, Hubbard L, Cruickshanks KJ, Palta M: 
The relation of retinal vessel caliber to the incidence and progression of 
diabetic retinopathy: XIX: the Wisconsin Epidemiologic Study of Diabetic 
Retinopathy. Arch Ophthalmol 2004, 122(1)76-83. 

1 2. Kifley A, Wang JJ, Cugati S, Wong TY, Mitchell P: Retinal vascular caliber, 
diabetes, and retinopathy. Am J Ophthalmol 2007, 1 43(6):1 024-1 026. 

1 3. Nguyen TT, Wang JJ, Sharrett AR, Islam FM, Klein R, Klein BE, Cotch MF, 
Wong TY: Relationship of retinal vascular caliber with diabetes and 
retinopathy: the Multi-Ethnic Study of Atherosclerosis (MESA). Diabetes 
Care 2008, 31(3)544-549. 

14. Roy MS, Klein R, Janal MN: Retinal venular diameter as an early indicator 
of progression to proliferative diabetic retinopathy with and without 
high-risk characteristics in African Americans with type 1 diabetes mellitus. 
Arch Ophthalmol 201 1, 129(1):8-15. 

1 5. Murray CD: The Physiological Principle of Minimum Work: I. The Vascular 
System and the Cost of Blood Volume. Proc Natl Acad Sci USA] 926, 
12(3)207-214. 

16. Chapman N, DeH'omo G, Sartini MS, Witt N, Hughes A, Thorn S, Pedrinelli R: 
Peripheral vascular disease is associated with abnormal arteriolar 
diameter relationships at bifurcations in the human retina. Clin Sci (Lond) 
2002, 103(2):111-116. 

1 7. Stanton AV, Wasan B, Cerutti A, Ford S, Marsh R, Sever PP, Thorn SA, Hughes 
AD: Vascular network changes in the retina with age and hypertension. 
JHypertens 1995, 13(12 Pt 2):1 724-1 728. 



18. Witt N, Wong TY, Hughes AD, Chaturvedi N, Klein BE, Evans R, McNamara M, 
Thorn SA, Klein R: Abnormalities of retinal microvascular structure and 
risk of mortality from ischemic heart disease and stroke. Hypertension 
2006, 47(5)575-981. 

1 9. King LA, Stanton AV, Sever PS, Thorn SA, Hughes AD: Arteriolar length-diameter 
(L:D) ratio: a geometric parameter of the retinal vasculature diagnostic of 
hypertension. J Hum Hypertens 1996, 10(6)417-418. 

20. Chapman N, Mohamudally A, Cerutti A, Stanton A, Sayer AA, Cooper C, 
Barker D, Rauf A, Evans J, Wormald R, Sever P, Hughes A, Thorn S: Retinal 
vascular network architecture in low-birth-weight men. J Hypertens 1997, 
15(12 Pt 1):1449-1453. 

21 . Patton N, Pattie A, MacGillivray T, Aslam T, Dhillion B, Gow A, Starr JM, 
Whalley U, Deary IJ: The association between retinal vascular network 
geometry and cognitive ability in an elderly population. Invest 
Ophthalmol Vis Sci 2007, 48(5):1 995-2000. 

22. Cheung N, Donaghue KC, Liew G, Rogers SL, Wang JJ, Lim SW, Jenkins AJ, 
Hsu W, Li Lee M, Wong TY: Quantitative assessment of early diabetic 
retinopathy using fractal analysis. Diabetes Care 2009, 32(1 ):1 06-1 10. 

23. Sasongko MB, Wang JJ, Donaghue KC, Cheung N, Benitez-Aguirre P, Jenkins 
A, Hsu W, Lee ML, Wong TY: Alterations in retinal microvascular geometry 
in young type 1 diabetes. Diabetes Care 2010, 33(6):1 33 1 -1 336. 

24. Benitez-Aguirre P, Craig ME, Sasongko MB, Jenkins AJ, Wong TY, Wang JJ, 
Cheung N, Donaghue KC: Retinal vascular geometry predicts incident 
retinopathy in young people with type 1 diabetes: a prospective cohort 
study from adolescence. Diabetes Care 2011, 34(7):1 622-1 627. 

25. Aldington SJ, Kohner EM, Meuer S, Klein R, Sjolie AK: Methodology 
for retinal photography and assessment of diabetic retinopathy: 
the EURODIAB IDDM complications study. Diabetologia 1995, 
38(4):437-444. 

26. Al-Diri B, Hunter A, Steel D, Habib M: Manual measurement of retinal 
bifurcation features. Conf Proc IEEE Eng Med Biol Soc 2010, 2010:4760-4764. 

27. Zamir M: Nonsymmetrical bifurcations in arterial branching. J Gen Physio 
1978, 72(6):837-845. 

28. Rassam SM, Patel V, Brinchmann-Hansen O, Engvold O, Kohner EM: Accurate 
vessel width measurement from fundus photographs: a new concept. 

Br J Ophthalmol 1994, 78(1)24-29. 

29. Brinchmann-Hansen O, Engvold O: Microphotometry of the blood column 
and the light streak on retinal vessels in fundus photographs. Acta 
Ophthalmol Scand 1986, Supplement 179:9-19. 

30. Gregson PH, Shen Z, Scott RC, Kozousek V: Automated grading of venous 
beading. Comput Biomed Res 1995, 28(4)291-304. 

31. Zhou L, Rzeszotarski MS, Singerman U, Chokreff JM: The detection and 
quantification of retinopathy using digital angiograms. IEEE Trans Med 
Imaging 1994, 13(4):61 9-626. 

32. Lowell J, Hunter A, Steel D, Basu A, Ryder A, Kennedy RL: Measurement of 
retinal vessel widths from fundus images based on 2-D modeling. 

IEEE Trans Med Imaging 2004, 23(10):! 196-1204. 

33. Habib M: Assessment of Retinal Vascular Geometry in Normal and Diabetic 
subjects. MD thesis. University of Newcastle Upon Tyne, School of 
Neuroscience: 2012. 

34. Rogers SL, Tikellis G, Cheung N, Tapp R, Shaw J, Zimmet PZ, Mitchell P, 
Wang JJ, Wong TY: Retinal arteriolar caliber predicts incident retinopathy: 
the Australian Diabetes, Obesity and Lifestyle (AusDiab) study. Diabetes 
Care 2008, 31(4)761-763. 

35. Nguyen TT, Wang JJ, Islam FM, Mitchell P, Tapp RJ, Zimmet PZ, Simpson R, 
Shaw J, Wong TY: Retinal arteriolar narrowing predicts incidence of 
diabetes: the Australian Diabetes, Obesity and Lifestyle (AusDiab) Study. 
Diabetes 2008, 57(3)536-539. 

36. Fanucci E, Orlacchio A, Pocek M, Magrini A, Salomoni E: Optimal branching 
of human arterial bifurcations. Invest Radiol 1990, 25(l):62-66. 

37. Klein R, Klein BE, Moss SE, Wong TY, Hubbard L, Cruikshanks KJ, Palta M: 
Retinal vascular abnormalities in persons with type 1 diabetes: the 
Wisconsin Epidemiologic Study of Diabetic Retinopathy: XVIII. 
Ophthalmology 2003, 1 1 0(1 1 ):21 1 8-21 25. 

38. Leung H, Wang JJ, Rochtchina E, Wong TY, Klein R, Mitchell P: 
Dyslipidaemia and microvascular disease in the retina. Eye (Lond) 2005, 
19(8):861-868. 

39. Bronson-Castain KW, Bearse MA Jr, Neuville J, Jonasdottir S, King-Hooper B, 
Barez S, Schneck ME, Adams AJ: Early neural and vascular changes 

in the adolescent type 1 and type 2 diabetic retina. Retina 2012, 
32(0:92-102. 



Habib ef al. BMC Ophthalmology 2014, 14:89 
http://www.biomedcentral.com/1471-2415/14/89 



Page 11 of 1 1 



40. Hubbard LD, Brothers RJ, King WN, Clegg LX, Klein R, Cooper LS, Sharrett 
AR, Davis MD, Cai J: Methods for evaluation of retinal microvascular 
abnormalities associated with hypertension/sclerosis in the 
Atherosclerosis Risk in Communities Study. Ophthalmology 1999, 
106(12):2269-2280. 

41. Porta M, Allione A: Current approaches and perspectives in the medical 
treatment of diabetic retinopathy. Pharmacol Ther 2004, 103(2):167-177. 

■ 

doi:10.1 186/1471-2415-14-89 

Cite this article as: Habib ef at The association between retinal vascular 
geometry changes and diabetic retinopathy and their role in prediction 
of progression - an exploratory study. BMC Ophthalmology 2014 14:89. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at /^'\ n! _, ul _-| r Q r,tr=l 

www.biomedcentral.com/submit ammBa central 



V 



